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FOREWORD 


The  research  reported  here  was  sponsored  by  the  Arnold  Engi¬ 
neering  Development  Center,  Air  Force  Systems  Command,  Arnold 
Air  Force  Station,  Tennessee,  under  Program  Element  61102F, 

Project  8951,  Task  06.  Air  Force  program  No.  61102F  is  entitled 
Research  and  Experimental  Study  of  Surface  Catalysis  in  Nonequi¬ 
librium  Flows  which  is  devoted  to  developing  a  probe  with  a  discon¬ 
tinuous  catalytic  surface  as  a  diagnostic  device  for  high  temperature 
wind  tunnels. 

The  research  was  performed  in  the  Aerodynamic  Research  De¬ 
partment  of  the  Cornell  Aeronautical  Laboratory,  Inc.  at  Buffalo, 

New  York,  under  Contract  No.  F40600-69-C-0005  during  the  period 
September  1967  to  September  1969  under  the  technical  cognizance  of 
Lt.  Vincent  Rocco,  AELR.  This  particular  theoretical  study  was 
initiated  to  provide  a  basis  for  determining  the  onset  of  gas-phase 
reactions  in  the  boundary  layer  on  the  probe.  Hence,  the  theory  pro¬ 
vides  a  basis  for  defining  the  range  of  applicability  for  a  catalytic 
probe.  This  report  is  also  identified  as  CAL  Report  No.  AF-2753-A-1 
and  was  released  for  publication  on  14  November  1969, 

This  technical  report  has  been  reviewed  and  is  approved. 

Vincent  A.  Rocco  Harry  L.  Maynard 

1st  Lt,  USAF  Colonel,  USAF 

Research  Division  Director  of  Plans 

Directorate  of  Plans  and  Technology 

and  Technology 
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ABSTRACT 

An  existing  theory  for  predicting  the  onset  of  dissociation  in  a  laminar 
boundary  layer  is  generalized  and  extended  to  account  for  dissociation  in  the 
boundary  layer  and  gas-phase  recombination  near  a  cold  wall.  The  existing 
theory  shows  that  dissociation  first  occurs  in  a  thin  layer  centered  on  the 
maximum  temperature  in  the  boundary  layer.  The  present  extension  shows 
that  for  cold  walls,  there  will  be  a  dense  layer  near  the  wall  with  relatively 
long  particle  residence  times  to  promote  recombination.  Hence,  at  least 
initially,  the  reacting  boundary  layer  on  a  cold  wall  is  characterized  by  a 
dissociation  layer  embedded  in  the  boundary  layer  and  a  recombination  layer 
near  the  wall.  The  theory  is  derived  without  specifying  the  gas  model,  and 
numerical  results  are  given  which  are  based  upon  a  realistic  gas  model.  It 
is  shown  that  the  effects  of  dissociation  and  recombination  tend  to  compete, 
and  under  certain  conditions  the  effects  at  the  wall  tend  to  cancel. 
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I.  INTRODUCTION 

The  occurrence  and  implications  of  thermochemical  nonequilibrium  in 

hypersonic  flowfields  has  been  widely  recognized,  and  there  is  a  large  body 

of  theoretical  and  experimental  literature  addressed  to  certain  aspects  of  the 

problem.  The  major  emphasis  has  been  on  chemical  relaxation  in  the  inviscid 

1-3 

flowfield,  and  theoretical  methods  are  available  for  idealized  gas  flows  , 

4-7 

for  comprehensive  models  of  reacting  air  flows  ,  and  for  scaling  the  initial 

8  9 

relaxation  processes  '  .  The  importance  of  thermochemical  relaxation  in 
viscous  flowfields  has  been  acknowledged,  but  the  body  of  literature  dealing 
with  this  aspect  is  relatively  sparse.  Most  of  the  literature  on  this  topic 
reports  machine  solutions  for  specific  isolated  problems^’  *  ,  and  while 
these  offer  important  information  on  reacting  viscous  layers,  it  is  difficult 
to  scale  these  results  to  other  problems  of  practical  interest. 

One  problem  not  treated  fully  in  the  literature  is  chemical  relaxation 
in  the  boundary  layer  of  a  slender  body.  This  problem  is  of  potential  impor¬ 
tance  in  hypersonic  wind  tunnel  experiments  with  a  cold  model  in  which  the 
ambient  flow  is  not  equilibrated,  and  consequently  contains  a  significant  con¬ 
centration  of  atomic  species.  Chemical  relaxation  could  occur  in  two  ways. 

As  noted  by  Rae^,  viscous  heating  will  produce  temperatures  in  the  boundary 
layer  that  exceed  those  in  the  inviscid  flowfield.  Consequently,  dissociation 
would  occur  initially  in  a  thin  layer  centered  on  the  maximum  temperature  in 
the  boundary  layer.  Rae's  analysis  can  be  applied  to  predict  the  onset  and 
consequences  of  this  relaxation.  For  a  cold-wall  model,  there  is  a  second 
and  competing  recombination  effect  which  could  take  place.  In  the  vicinity  of 
the  wall,  the  low  gas  temperature  coupled  with  the  low  gas  velocities  could 
produce  the  high  densities  and  long  particle  residence  times  required  for 
recombination.  Consequently,  recombination  could  occur  initially  in  a  thin 
layer  near  the  wall.  These  two  relaxation  processes  could  begin  sequentially 
or  simultaneously,  depending  on  the  ambient  conditions,  thereby  altering  the 
energy  transport  to  the  surface.  These  effects  could  be  important  in  routine 
experiments  to  determine  heat  transfer  rates  for  particular  configurations. 

It  would  be  of  crucial  importance  in  applications  of  the  slender  catalytic 
probe  being  developed  in  this  contract. 
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The  purpose  of  this  report  is  to  present  an  approximate  theoretical 
method  for  predicting  the  onset  of  chemical  relaxation  processes  in  the  laminar 
boundary  layer  on  a  two-dimensional  flat  plate  or  wedge.  The  theory  is  an  ex¬ 
tension  of  Rae's  theory  to  include  recombination  near  the  wall  and  a  generali¬ 
zation  in  the  sense  that  Liighthill's  idealized  gas  model  is  not  used. 

Rae's  analysis  is  paralleled  in  the  next  section  to  arrive  at  the  most 
general  theoretical  result  without  specifying  a  gas  model.  Following  that> 
consideration  is  given  to  the  most  important  recombination  and  dissociation 
reactions  in  high  temperature  air  to  arrive  at  numerical  results. 
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II.  THEORY 


The  theoretical  development  begins  by  writing  down  the  boundary  layer 
equations  for  a  reacting  gas  and  transforming  these  using  the  Levy-Lees  co- 

13 

ordinates.  The  boundary  layer  equations  are  available  in  any  standard  text  ’ 
and  are  not  repeated  here.  It  is  assumed  that  both  the  Prandtl  number  and  the 
Schmidt  number  are  unity  and  that  the  viscosity  of  the  gas  varies  linearly  with 
temperature.  Moreover,  it  is  assumed  that  variations  in  the  specific  heat 
through  the  boundary  layer  are  negligibly  small,  and  that  the  difference  in  the 
enthalpy  of  the  atomic  and  molecular  species  is  small  in  comparison  with  the 
heat  of  formation.  These  approximations  are  consistent  with  the  spirit  of  the 
analysis;  to  identify  the  initial  relaxation  processes. 


Using  the  transformation 

5  „n„u.Jz  1  ■  Jkffdy  (1) 

the  boundary  layer  equations  become 


t  ffn  - 

z  ■  -Pf  f ti ) 

(2) 

wr,  * 

f  -  £  T,) 

Mw  C+  U* 

(3) 

+  + 

ZS  uT, 

=  25  K  »r  " 

P-W  9 

(4) 

If  it  is  assumed  that  boundary  layer  displacement  effects  are  negligible,  then 

terms  involving  fg  can  be  neglected.  In  this  case,  Eq.  (2)  reduces  to  the 

Blasius  equation,  and  the  solution  can  be  taken  from  any  standard  text,  for 
* 

example  Ref.  15.  The  energy  and  species  equations  are  made  nondimen sional 
by  dividing  by  a  reference  temperature,  lAe* / Cj,  ,  and  the  atom  mass  fraction 
at  the  edge  of  the  boundary  layer,  .  An  additional  transformation  is  made, 

c  .  — S-  - 

?w  u*  T 

*It  should  be  noted  that  Ref.  15  provides  a  solution  for  the  equation, 

and  if  these  results  are  used,  the  appropriate  numerical 
transformation  must  be  made. 
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where  T  is  a  characteristic  time  to  be  determined  with  the  chemical  rate 
equations.  With  these  manipulations,  the  energy  and  the  species  equations 
become 


a  0L.iL  uf'  X 

s„  ♦  / s,  -  2ff,sr  -c  •  -  2  nf  75: 


A„  +  n,  -  2f  t,i. 


hfiT 

po£e 


(5) 

(6) 


The  method  used  by  Rae  to  obtain  solution  to  Eq.  (5)  and  (6)  is  to 
assume  that  the  dimensionless  temperature  and  atpm  fraction  can  be  expressed 
in  a  series  form 


G(r)£)  -  +  C  +  "• 


and  to  substitute  these  in  Eq.  (5)  and  (6)  to  obtain  a  system  of  ordinary 
differential  equations.  Using  just  the  first  two  terms  in  each  series  expan¬ 
sion,  the  resulting  system  of  equations  is 


®.„  -  V  ■  * 

with  the  boundary  conditions, 

*]  -*  0  :  0o  — • -  }  A,—*  1.0 

n  —  00  :  $  — -  ,  A„  — «■  1.0 


(7) 

(8) 


(9) 

(10) 
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0  —  0  •'  —  fl  ,  A,,  —  o 
f\  — ►  oo  :  Gi  — *  0  t  A,  — ►  0 

Equations  (7)  and  (8)  are  the  governing  relations  for  frozen  flow  and  their 
solutions  are 

jf,  (1  +2<9«) 

A.  =  i.o 


(ID 

(12) 


Equations  (9)  and  (10)  govern  the  initial  departures  from  the  frozen 
flow,  and  because  of  the  series  form  assumed  for  the  solution,  the  production 
term,  - ,  is  evaluated  at  the  frozen  conditions.  This  single  feature 
makes  an  analytic  solution  possible.  Continuing  to  follow  Rae's  analysis, 
the  method  of  solution  is  to  make  the  standard  substitution, 

0f  -  Slrj)  U(q)  A,  *  SO])  VJ(^) 

to  arrive  at  the  following  system  of  equations. 


+  $ 

- 

2  f,S  =  0 

(13) 

4*  *  2  °n 

5, 

5 

L  +  *  f 

u*  ( 

e 

e<*c  ° 

(14) 

5, 

S 

1  ♦  fv,-- f 

/ 

). 

(15) 

with  the  boundary  conditions 

jfj  s  0  S  (o)  *  1  j  5^  (o)  =  0,  0(0)  *  0  ,  3  0 

iq  — •>  co  ;  U  =  W  -  0 

The  boundary  conditions  on  S(*j)  are  arbitrary  and  were  selected  by  Rae 
for  convenience.  Those  on  WOj)  and  l/(^)  are  chosen  to  be  consistent  with 
the  boundary  conditions  for  Eq.  (9)  and  (10).  Rae  has  solved  Eq.  (13)  numer¬ 
ically  and  tabulates  that  solution  in  Ref.  12  along  with  the  associated 
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function, 


f  ‘ct? 
s*<o 


cL  * 


where  the  prime  notation  indicates  differentiation  with  respect  to  the  dummy 
variable,  t  .  They  are  reproduced  here  in  Table  I  for  convenience.  The 
solutions  to  Eq.  (14)  and  (15)  are  obtained  in  a  conventional  manner  and  are 


U(rj) 


C  -  iTty)  G0)> 


W(q) 


C  - 


VO])  G(f}) 


v(t)  s (t)  r  u/.r 

nt )  lea« 


C  »  2  V(oo) 

Grty)  '  2. 


n>- 

vlt)l 

Set) 

'  uxLZ' 

l.  r 

VO®)] 

■ret) 

. 

S(t) 

r^r 

oL 

f'U) 

1 

0 

dL  t 


6  =  +•  5  S(^|)  L/O]) 

=  r  *  C  soj)  wty) 


rit)  s it)  \  or- Z]  ,  ,  1 

r«>  U].dtJ 

(16) 

oi  t 

(17) 

0 

(18) 

(19) 

(20) 

(21) 

The  above  relations,  combined  with  Eq.  (11)  and  (12)  provide  a  complete 
solution  to  the  problem. 


Two  quantities  of  interest  are  the  atom  concentration  at  a  noncatalytic 
wall  and  the  convective  heat  transfer  to  a  noncatalytic  wall.  Noting  that  for  a 


transfer  are 


K  - 

* 

Uer  *  toe  1 

-  1  +  c  - 

ot* 

x> 

$ 

-  i  -  *-5 

V 

irioo ) 


u*T  l/£  +  0, 


(22) 


(23) 
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where  denotes  the  convective  heat  transfer  when  there  are  no  chemical 
reactions  present.  It  can  be  seen  that  the  solution  depends  upon  evaluating  the 
constant*  C  ,  which  in  turn  requires  integrating  the  production  term  across 
the  boundary  layer.  A  similar  integration  involving  the  production  term  is 
also  required  to  obtain  the  species  and  energy  profiles. 


A  general  rate  equation  for  dissociation  and  recombination  is 


(24) 


where  the  subscripts,  m.  and  A  ,  refer  to  the  molecular  and  atomic  state 
of  the  reacting  species,  and  the  subscript,  *  ,  refers  to  the  inert  species 

that  enter  as  collision  partners.  The  general  relations  for  the  equilibrium 
constants  and  the  forward  and  reverse  reaction  rates  are 


i y  s'  T  " n<:  „  "  "rpl'r  JL  S'  T  Xp/t  A 

Kt-  CaT  e  -kfti  -  C4i  T  e  itr,i 


K, 


-Cr\T 


-n 


ri 


and  where  Xp  is  the  characteristic  temperature  for  dissociation. 


Restricting  attention  to  the  case  Tw  ^  Te  ,  it  can  be  seen  that 
integration  of  Eq.  (24)  will  produce  two  parts.  In  the  vicinity  of  the  wall,  the 
temperature  will  be  low  relative  to  the  rest  of  the  boundary  layer  and  the 
dissociation  component  of  Eq.  (24)  will  be  exponentially  small  in  comparison 
with  the  recombination  component.  Since  the  present  analysis  is  concerned 
only  with  the  initial  relaxation  process,  it  can  be  argued  further  that  the  re¬ 
combination  term  will  be  important  only  in  the  coldest  part  of  the  boundary 
layer;  a  narrow  region  immediately  at  the  wall.  The  same  arguments  apply 
to  the  dissociation  component,  as  pointed  out  by  Rae.  Initially  at  least, 
dissociation  will  be  important  only  in  a  narrow  region  centered  on  the  maxi¬ 
mum  temperature  in  the  boundary  layer,  and  the  recombination  term  can  be 
neglected  in  effecting  this  part  of  the  integration.  As  a  consequence,  Eq.  (18) 
takes  on  the  form,  C  =  C*  +  Cv  ,  where  CR  and  C,  are  the  recombination 
and  the  dissociation  components  of  the  integration  and  are  given  by 


c* 

C, 


-  4  V(oo)  X 

A-  v(e»)T 


1. 1 

ft 


I  - 

7  - 


Wt)l  S(t;  ,  ?«■„  m,  j*; 

V0*»J  fit)  K  J  i  mi  CL*  r>i 

VL t?  ]  _Sct)  P_ 

fit)  7  wTi  "otl  ° 


(18  a) 
(18b) 
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The  evaluation  of  Eq.  (18)  can  proceed  in  two  separate  steps,  first  considering 
the  recombination  reaction  near  the  wall  and  then  the  dissociation  reaction  in 
the  boundary  layer. 


WALL  RECOMBINATION 


The  method  used  to  evaluate  Eq.  (18a)  is  based  on  the  fact  that 
attention  is  restricted  to  a  thin  region  near  the  wail.  In  this  region,  the 
following  approximations  are  valid. 

^  '(*7>  »  f  1  ,  <50(r|)  ®  6W  +  [l  4  Z  Ot  -Z  flj  f%), ;  ,  virj)  «  ^  ,  S  =  1.0 


Using  these  and  the  general  expression  for  the  reverse  reaction  rate,  Eq.  (18) 
becomes  after  some  algebraic  manipulation, 


4-X  V(eo)  f 

2 

Vt„1 

r  ^ 
L  m; 

[r«>]s  1 

i  _  i» 

7^  T*.  - 

h-2>.  A, 

and 

A 

r,i  (T*,)  indicates 

CL* 


(1  +  Vl-i) 


evaluated  at  the  wall  temperature.  The  above  expression  neglects  (6>e) 
in  comparison  with  U9„)  **  ,  which  is  consistent  with  the  cold  wall 

approximation  introduced  earlier. 


(25) 


DISSOCIATION  REACTION 

A  comprehensive  discussion  of  the  technique  for  integrating  Eq.  (18) 

for  the  dissociation  reaction  is  given  in  Ref.  12.  Briefly,  it  is  based  upon 

the  previous  observation  that  the  reaction  is  confined  to  a  thin  region  centered 

on  the  maximum  temperature  in  the  boundary  layer.  This  makes  it  possible 

1Z  16 

to  approximate  the  integral  using  the  Laplace  method  ’  .  This  method 

states  that  if  in  the  integral 

/  -  8 

\/Ct)  e  cLt 

Vet)  is  bounded  in  the  interval  from  a-  to  t  ,  and  if  there  exists  a  point 
in  this  interval  where  oJ1  (tA)  *  0  ,  u >  0  ,  then  for  B  >  0  ,  the  integral 
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can  be  approximated  by 


2  t r 


8  or  "(t*) 


vtt.)  e 


-  8nf  (t  *  ) 


Applying  this  to  Eq.  (18b),  and  after  some  algebraic  manipulation,  the  rela¬ 
tion  for  the  dissociation  component  of  the  constant  becomes 

vn.t  of.*  (26) 


4/fl  vCco)  T  \/z  Tf  %  7  SLL  -fc.(T) 

V  e  rD  4*  vn.:  cl.  ^ 


where  the  notation,  t  (T,)  indicates  the  forward  reaction  rate  is 
evaluated  at  the  temperature,  Tm  , 

<5*  "  *  (j  +  *  *,>' 

As  indicated  in  Ref.  12,  the  above  method  applies  only  for  6e  ~  <  £ 

The  constant  tR  in  Eq.  (26)  is  defined  as 


(27) 


„  „  s<-^  rf .  ZHt) 

i-  vloo> 

It  can  be  taken  to  be  unity  with  errors  of  less  than  15%  for  0.  18  <  (,Oe  -  0J<-  0.46 
which  covers  much  of  the  range  of  interest. 


The  parameter  governing  the  species  concentrations  at  the  wall  and 
the  convective  heat  transfer  to  the  wall  is  f  CR*  Cft 1  .  This  follows  from 

Ue?  I  V(eo)  J 

Eq.  (25)  and  (26) . 


Z  C 
Ue  T  vloo) 


(28) 


The  two  terms  in  the  braces  can  be  identified  respectively  as  •—  and 

Tp 

where  L  is  a  characteristic  relaxation  time  and  the  subscripts  D  and  R 
refer  to  dissociation  and  recombination.  The  competing  effects  of  gas-phase 
dissociation  and  recombination  are  apparent.  It  is  possible,  at  least  in 
principle,  to  have  significant  dissociation  and  recombination  occurring  within 
the  boundary  layer  without  any  observable  effect  at  the  surface.  This  would 
correspond  to  the  two  terms  within  the  braces  being  equal  and  non-zero. 
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GAS  MODEL 

A  simplified  gas  model  for  air  will  be  used  to  assess  the  importance 
of  gas-phase  dissociation  and  recombination  in  the  boundary  layer  of  a  flat- 
plate.  The  primary  interest  here  is  on  experiments  in  a  hypersonic  wind 
tunnel  which  is  characterized  by  a  nonequilibrium  flow  in  the  nozzle  expansion. 
Consequently,  the  ambient  flow  would  be  a  mixture  of  atomic  and  molecular 
oxygen,  nitric  oxide,  and  molecular  nitrogen.  Atomic  nitrogen  will  not  be 
considered  here  because  it  is  energetically  important  only  at  very  high  tem¬ 
peratures.  The  two  reactions  considered  are 

Ax 

(a)  02  +  yn  — ~  0  +  0  +  yn 


(b) 


blO^  +  Yn. 


NO  +  0  +  Yn. 


The  first  is  the  dominant  reaction  for  producing  oxygen  atoms,  but  the  second 
is  an  important  reaction  for  oxygen  atom  consumption  at  low  temperatures. 
Rigorously,  a  third  reaction  should  be  included,  NO  +  0 &  — £»AI0l+-  0  ,  because 
the  reverse  reaction  can  also  be  important  at  low  temperatures.  However, 
only  the  species  present  in  the  ambient  flow  enter  into  the  results,  and  the 
only  source  of  N0:  is  the  products  of  the  (b)  reaction.  Hence  the  above  re¬ 
action  must  be  neglected  in  the  present  analysis.  The  rate  equation  for  the 
two  reactions  listed  is 


eta. 

alt 


YTU. 


Z,  yn.  /  V  \ 

Zkf''k  ^ 

i  } 

p  *<*.**.  A 

>n.  m*. 


where  the  subscripts,  0  ,  ^  ,  and  NO  refer  to  that  species,  and 

the  subscript  >n.  refers  to  the  inert  third  body. 
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The  characteristic  relaxation  time  identified  in  Eq.  (28)  can  be  cast 
into  the  following  form. 

*  „  ,a  V r»  r  ^  a  n  ) 


(29  b) 


'  -  IQ  V-  r  a-:  J.  fT^ 

*7  ^  *  -  -  f*:  4  (29a) 

J8  12  \*  V  -4r.;CTw) 

_J_  _  I  m.  y  U  m-i  a.  + 1) 

T«  “  Ltfl.  £-)  [i.|  **]•  <«b 

The  reaction  rates  were  evaluated  using  the  data  given  in  Refs.  17-20  and 
summarized  in  Table  II.  The  resulting  numerical  expressions  for  use  in 
Eq.  (29)  above  are  _ 


where  1K  =  300  °K  is  an  arbitrary  reference  temperature.  The  character¬ 
istic  relaxation  times  have  been  numerically  evaluated  assuming  the  bracketed 
terms  in  the  above  reaction  rates  are  unity.  These  are  plotted  in  Figs.  1 
and  2.  It  should  be  noted  the  relaxation  time  in  Figs.  1  and  2  is  overestimated 
because  of  the  above  assumption.  For  typical  conditions,  the  error  amounts 
to  10%  to  25%  and  it  can  be  corrected,  if  desired,  by  multiplying  the  ordinate 
by  the  appropriate  value  of  the  bracketed  quantity. 

It  will  be  noted  that  Fig.  1  was  evaluated  for  ^/T,  =  0.  The  wall 

temperature  ratio  enters  into  the  determination  for  the  temperature  govern¬ 
ing  the  dissociation,  7*  ,  and  it  is  a  relatively  weak  function  of  the  wall 

temperature.  The  wall  temperature  ratio  was  taken  to  be  zero  for  this  reason 
and  because  present  interests  are  primarily  on  the  cold  wall  case.  It  will  be 
noted  in  Fig.  1  that  the  characteristic  time  for  dissociation  goes  through  a 
minimum.  This  behavior  arises  because  of  the  pre- exponential  temperature 
dependence  in  the  rate  equation.  It  is  of  academic  interest  only  because  it 
corresponds  to  a  stagnation  temperature  of  T0  ^  10^  °K. 
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The  characteristic  time  for  recombination!  Fig.  2,  is  a  simple 
function  of  the  wall  temperature  ratio  and  the  Mach  number  at  the  edge  of  the 
boundary  layer.  It  was  evaluated  for  the  case  of  T^.  =  300  °K,  but  the  results 

can  be  corrected  to  any  other  wall  temperature  by  dividing  the  ordinate  by 

(  T„  )3'2. 

It  was  noted  earlier  that  it  would  be  possible  to  have  signification 
dissociation  and  recombination  occurring  in  the  boundary  layer,  but  because 
the  effects  tend  to  compete,  there  might  be  no  net  effect  at  the  wall.  The 
following  conditions  illustrate  this  effect.  For  ^  =  &-ox  “ 

and  Me  *  4  ,  the  characteristic  times  for  dissociation  and  recombination 

are  equal  for  Tp/T0  *  4  ,  and  are  zjr  -  i-  ^  *  1 03/see..  Using  a 

rj>  4 

typical  length  of  one  foot  and  a  typical  velocity  of  10  ft/  sec,  Eq.  (22)  shows 
that  each  effect  could  alter  the  atom  fraction  at  the  wall  by  about  10%.  However, 
because  the  characteristic  times  are  equal,  the  two  effects  just  cancel  and 
leave  the  wall  concentration  unaltered. 
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III.  CONCLUDING  REMARKS 

An  existing  theory  for  predicting  the  onset  of  dissociation  in  a 

12 

laminar  flat  plate  boundary  layer  has  been  extended  to  account  for  gas- 
phase  recombination  near  the  wall.  The  basis  for  this  extension  is  the 
observation  that  near  a  cold  wall,  the  combination  of  the  low  gas  tempera¬ 
ture  and  low  velocity  could  produce  the  high  density  and  relatively  long 

particle  residence  time  necessary  for  gas -phase  recombination.  In  affecting 

12 

this  extension,  Rae's  theory  is  generalized  to  allow  for  a  more  realistic 
gas  model  of  air.  The  extension  to  account  for  recombination  centers  on 
approximating  the  boundary  layer  profiles  near  the  wall  in  order  to  integrate 
the  species  production  term  in  the  region  where  recombination  is  important. 

Numerical  results  were  obtained  to  illustrate  the  effects  of  the  gas- 
phase  reactions.  These  results  show  that  dissociation  can  occur  in  the 
boundary  layer  while  gas-phase  recombination  is  occurring  near  the  wall. 
These  effects  tend  to  compete,  and  under  certain  conditions,  significant 
dissociation  and  recombination  can  occur  with  no  net  change  in  atom  con¬ 
centration  or  heat  transfer  to  the  wall. 
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TABLE  I 

FUNCTIONS  IN  THE  BOUNDARY  LAYER  SOLUTION 
(FROM  REF.  12) 


1 

vt» p 

f  V)> 

0. 

1 .000000 

0. 

0.46960 

0.  1 

1.000157 

0.04696 

0.46956 

0.2 

1.001252 

0.09385 

0.46931 

0.3 

1.004227 

0. 14051 

0.46861 

0.4 

1.010021 

0. 18667 

0.46725 

0.5 

1.019576 

0.23198 

0.46503 

0.6 

1.033840 

0.27599 

0.46173 

0.7 

1.053762 

0.31822 

0.45718 

0.8 

1.080302 

0.35817 

0.45119 

0.9 

1.  114429 

0.39539 

0.44303 

1.0 

1 . 157121 

0.42950 

0.43438 

1.  1 

1.209372 

0.46021 

0.42337 

1.2 

1.272191 

0.48736 

0.41057 

1.3 

1.346600 

0.51096 

0.39599 

1.4 

1.433639 

0.53110 

0.37969 

1 .5 

1  .534358 

0.54799 

0.36181 

1 .6 

1 .649818 

0.56194 

0.34249 

1 .7 

1.781086 

0.57327 

0.32195 

1.8 

1.929224 

0.58236 

0.30045 

1.9 

2.095288 

0.58953 

0.27825 

2.0 

2.280312 

0.595  14 

0.25567 

2.  1 

2.485301 

0.59946 

0.23301 

2.2 

2.711216 

0.60276 

0.21058 

2.3 

2.958970 

0.60526 

0. 18867 

2.4 

3.229407 

0.60712 

0.  16756 

2.5 

3.523301 

0.60851 

0.  14748 

2.6 

3.841343 

0.60953 

0.12861 

2.7 

4.184138 

0.61028 

0.11112 

2.8 

4.552198 

0.61082 

0.09511 

2.9 

4.945944 

0.61121 

0.08064 

3.0 

5.365708 

0.61149 

0.06771 

3. 1 

5.811736 

0.61169 

0.05631 

3.2 

6.284194 

0.61183 

0.04637 

3.3 

6.783178 

0.61193 

0.03782 

3.4 

7.308724 

0.61200 

0. 03054 

3.5 

7.860816 

0.61205 

0.02442 

3.6 

8.439400 

0.61208 

0.01933 

3.7 

9.044389 

0.61211 

0.01515 

3.8 

9.175679 

0.61212 

0.01176 

3.9 

10.333152 

0.61213 

0.00904 

4.0 

11.016685 

0.61214 

0.00687 

4.1 

11.726157 

0.61214 

0.00518 

4.2 

12.461451 

0.61215 

0.00386 

4.3 

13.222460 

0.61215 

0.00285 

4.4 

14.009087 

0.61215 

0.00208 

4.5 

14.821246 

0.61215 

0.00151 

4 . 6 

15.658863 

0.61215 

0.00108 

4.7 

16.521875 

0.61215 

0.00077 

4.8 

17.410231 

0.61215 

0.00054 

4.9 

18.323887 

0.61215 

0.00037 

5.0 

19.262809 

0.61215 

0.00026 

5.  1 

20.226971 

0.61215 

0.00018 

5.2 

21.216352 

0.61215 

0.00012 

5.3 

22.230935 

0.61215 

0.00008 

5.4 

23.270708 

0.61215 

0.00005 

5.5 

24.335662 

0.61215 

0.00003 

5.6 

25  .425791 

0.61215 

0.00002 

5.7 

26.541089 

0.61215 

■  4  ■ 
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TABLE  II 


RATE  CONSTANTS  FOR  THE  GAS  MODEL  OF  AIR 


*,  »  C,T-n‘  e'T*'T 


Cr  T’*' 


0Z  +  M 


0  +  0  +  Y/L 


(b)  NO  +  M-< — NO  *  0  +  nt 


Reaction 

Cf 

— 

c  k 

Ref. 

a 

O 

6.4  x  1023 

2.0 

20 

5.3  x  10 

3/2 

17 

°2 

1.9  x  1021 

3/2 

1  .6  x  1018 

1.0‘ 

18 

N2 

4.8  x  1020 

3/2 

4  x  1017 

1.0 

19 

b 

All 

22 

5.6  x  10 

5/2 

20 

17 


(i/sec) 


Figure  1 
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AED 


Figure  2  THE  CHARACTERISTIC  RELAXATION  TIME  FOR  OXYGEN  RECOMBINATION 

IN  A  LAMINAR  BOUNOARY  LAYER 


/=  1  A  T  =  30Q°K 
hr 
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